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Scaling laws and universality are often associated with systems exhibiting emergent phenomena 
possessing a characteristic energy scale. We report nonequilibrium transport measurements on two 
different types of single-molecule transistor (SMT) devices in the Kondo regime. The conductance 
at low bias and temperature adheres to a scaling function characterized by two parameters. This 
result, analogous to that reported recently in semiconductor dots with Kondo temperatures two 
orders of magnitude lower, demonstrates the universality of this scaling form. We compare the 
extracted values of the scaling coefficients to previous experimental and theoretical results. 



Scaling theories, most notably of the power law form, 
characterize a diverse array of naturally occurring phe- 
nomena. Seemingly unrelated systems with widely differ- 
ing microscopic details can be described by relations with 
identical scaling exponents when the underlying dynam- 
ics are similar. This property, known as universality, is 
relevant when considering systems driven out of equilib- 
rium by some perturbation. In a system with an emer- 
gent characteristic energy scale, examining the system 
response as a function of the perturbation scaled relative 
to that energy can reveal such universality 

The many-body Kondo state in quantum dots, as ob- 
served in a variety of microscopic implementations [l|, 0, 
S 0j 01, is an example of such a system. In these 
devices itinerant electrons in source and drain leads are 
coupled via tunneling barriers to a single magnetic im- 
purity. At zero source-drain bias {V), below an emergent 
Kondo energy scale, k^T^, the local moment of the im- 
purity is screened by the conduction electrons. As a re- 
sult of the screening process the (zero-bias) conductance, 
G, is enhanced at low temperatures (T < Tk). It has 
been documented that G(T, V ~ 0) in the equilibrium 
Kondo regime is well described by a universal tempera- 
ture dependence. 0, H, Q A nonzero applied bias {\V\ > 
0) drives the system into the nonequilibrium Kondo 



regime. 



121 G(T, V) exhibits a resonance peak cen- 
0. Theoretical treatments of idealized quan- 



tered at V 

tum dots have argued that G{T, V) in the single channel 
Kondo state may be described by an analogous function 
in bias voltage and temperature. However, there is dis- 
cussion regarding the energy scales and order to which 
this scaling will hold as well as the number of system 
specific coefficients required, their expecte d values, and 



the universality of these coefficients. Il3lll4l .ll5l. ll6l . ll7l .ll8j 
Recent experiments by Grobis et aL[l9| have shown that 
G(T, V) measured in a single channel GaAs quantum dot 
(Tk ~ 0.3 K) in the nonequilibrium regime is well de- 
scribed by a universal scaling function with two scaling 
parameters. 



In this Letter we test the universality of these re- 
sults, applying the analogous analytical approach to 29 
molecule-based devices with Kondo temperatures rang- 
ing from 35 K to 155 K. We find that the conductances 
of SMTs containing either Cgo @, 0, 21 1 or a transition 
metal complex are accurately described by the same 
scaled parameters as the GaAs dot[l^. We confirm the 
quadratic voltage and temperature dependence of G in 
the low energy limit. The values of the extracted scaling 
coefficients are quite consistent throughout the ensemble. 
We discuss possible explanations for the systematic dif- 
ferences between our coefficients and those inferred from 
previous experiments [l9l| and theoretical model predic- 
tions. 

Sample fabrication begins with e-beam lithographic 
patterning of gold nanoconstrictions with minimum 
widths between 90 and 120 nm, in arrays on an n"*" Si sub- 
strate with a 200 nm Si02 insulating layer. Immediately 
prior to being placed in a cryostat, the chip is cleaned and 
then treated with the molecules that are to be utilized 
as the active elements, as described previously @, [2^ . 
Each device in an array can be contacted individually 
using an Attocube Systems cryogenic probe station with 
a base temperature ~ 1.8 K. Electromigration in a cryo- 
genic environment is used to create a small gap in the 
nanowire, ideally allowing a molecule to bridge the bro- 
ken ends of the wire that will subsequently function as 
source and drain electrodes (Fig[T^ and fib) . [isj Geo and 
a Gu-containing transition metal complex (bis(2,5-di- 
[2]pyridyl-3,4-dithiocyanto-pyrrolate)Gu(II)) ( "complex 
1" ) , were independently utilized as the active elements 
of our SMTs. Both CgoBSIllI and complex l(2i,[24| 
have been successfully used as active elements in past 
investigations of breakjunction devices. 

Differential conductance, dl/dV, was measured as a 
function of V using quasi-DG lock-in techniques with 
an rms excitation voltage of 0.75 mV at 19.149 Hz. 
Post-electromigration transport properties were consis- 
tent with previous reported frequencies of occurrence. [2^ 



2 



A pronounced zero-bias conductance peak was observed 
in approximately 20% of devices, approximately a third 
of which were sufficiently stable at relevant biases and 
temperatures to enable a complete analysis. The proto- 
typical behavior of the single channel spin-^ Kondo effect 
is exhibited by a decreasing amplitude and broadening 
of the Kondo resonance as the temperature is increased. 
The temperature dependence of the amplitude has been 
well described by an empirical formQ derived from a fit 
to the renormalization group analysis. [Hj 

Gek(T, 0) = Go(l + (2i/-^ - 1){T/T^f)-' + Gb, (1) 

where s ~ 0.22 for a spin-i impurity. The total conduc- 
tance comprises the Kondo resonance and a background 
contribution, Gb, approximated as linear in applied bias, 
found via a fit to dl/dV vs V over a bias range that in- 
corporates the conductance minima outside the Kondo 
resonance. A constant is subtracted from the fit to en- 
sure that G(T, V^) - Gb > for all T and V. From Eq. 
Go is the height of the resonant peak, in the limit that 
T -> K: Go = G(0,0) - Gb. Similarly Tk is the 
temperature at which the Kondo resonance equilibrium 
conductance is equal to half its zero-temperature value: 
G{Tk, 0) — Gb = Go/2. From this description we extract 
values of the parameters Tk and Go (Fig. [If) . As T in- 
creases, deviations from Eq. ([!]) are expected as higher 
order transport processes become relevant. Our data is in 
good agreement with this form up to T/Tk = 0.4 — 0.75, 
depending on the device. Further references to G(T, V) 
denote conductance with the background component re- 
moved. 

Of the 29 devices analyzed over the full temperature 
range, 19 incorporated complex 1 and 10 utilized Ceo, 
with Tk ranging from 34 K to 155 K. While the Cgo- 
based devices had a lower Tk on average (TK,ave ~ 57 K 
for Ceo-based devices compared to TK,ave w 70 K for de- 
vices utilizing complex 1), they were typically less stable 
with increasing temperature and bias. The enhanced sta- 
bility of devices made with complex 1 is likely due to their 
thiocyanate end groups which can covalently bond to the 
gold electrodes. Irreversible changes in device configura- 
tion at higher temperatures limited our measurements to 
the region of T < Tk/2 for the majority of devices. 

If the Kondo resonance scales universally with respect 
to bias and temperature, then we can assume that the 
scaling function may be approximated by a series expan- 
sion. We test this idea by fitting the nonequilibrium con- 
ductance to the equation G(T, F) w Go — CT{k-s,T)^'^ — 
cv{eV)^^ . Py and Pt are the scaling exponents for bias 
and temperature, while cy and ct are expansion coef- 
ficients. We do not presume the theoretically predicted 
quadratic form of the scaling exponents and instead es- 
timate Py and Pt independently. To extract Py we fit 
G(T, 0) — G(T, V) to a power law form in the bias range 
\eV\ < fcBTK/2 for each trace in which T < Tk/5. To 
extract Pt this process is repeated for T < Tk/5 at fixed 
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FIG. 1: (color online), (a) Schematic of SMT measurement 
setup, (b) Transport path for conduction electrons in SMT 
device. Traces of dl /dV vs. V for several temperatures (c) 
from 1.8 K(black) to 53 K(maroon) for device A, using com- 
plex 1, and (d) from 1.8 K(black) to 46 K(maroon) for device 
B, using Ceo- Amplitude of Kondo resonance decreases with 
increasing temperature. Insets: Schematic of (c) complex 1 
and (d) Ceo- (e) dl/dV vs at T = 1.8 K for two sam- 
ples, demonstrating the disparate Kondo temperatures asso- 
ciated with different devices. Device C (blue) contains a Ceo 
molecule and D (red) contains complex 1. Note that the 
respective background contributions to the measured conduc- 
tance have not been subtracted from the plots in (c), (d), and 
(e). (f) Fit of equilibrium conductance versus temperature 
for devices A, B, C, and D, using equation [1] yielding Tk ~ 
52 K, 43 K, 35 K, and 105 K, respectively. 



values of V , again in the range \eV\ < /cbTk/2. Wc limit 
the T and V ranges because the nonequilibrium conduc- 
tance is expected to deviate from universal scaling as /cbT 
and eV approach /cbTk. Averaged over all devices, we 
find Pv = 1.97 ± 0.047 and Pt = 2.07 ± 0.18. Theo- 
retical treatments have found the leading order correc- 
tions to the conductance are quadratic in temperature 
and bias.[l^, [2^, [2^ The scaling exponents we find for 
both T and V are thus very near the predicted value of 
2. Henceforth we use Py and Pt equal to 2. The error 
bars for all extracted parameters in this study represent 
95% confidence intervals. 

Following the approach of Grobis et al. [l^ , we first ob- 
tain a scaling function of normalized temperature, T/T}^, 
and bias, eV/kBTK describing the nonequilibrium con- 
ductance at low energies. In the context of the Anderson 
impurity model [llll2a|. the Kondo resonance arises from 
the local density of states [131 of a magnetic impurity with 
charging energy Ec which is coupled to the electron bath 
of the leads; in a dot in the Kondo regime this is propor- 
tional to the conductance and may be expanded in the 
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FIG. 2: (color online). Extracted values of the scaling co- 
efficient (a) a and (b) 7 for each device. Open blue circles 
correspond to devices with complex 1. Filled red circles corre- 
spond to devices with Ceo. Dashed gray line indicates average 
value. Samples are numbered in order of increasing Tk. 



low energy limit as[l 



G(T,y) = Gek(T,0) 



1 - 



-(—] 



J 

(2) 

where a and 7 are the scaling coefficients we intend to 
extract. The constant ct ~ 4.92. This is set by the low 
temperature expansion of Eq. ([T|) in which Gek(2^, 0) « 
Gq{1 — ct{T /T^if')- The constants a, 7, and ct arc ar- 
ranged in this fashion so that Eq. ^ can be expanded at 
low temperature to a simplified power law form. Here, a 
and 7 are defined such that they are independent of the 
definition of Tk. Eq. ^ in the small T limit is of the 
formfli, [3 



G{T,0) -G{T,V) 
ctGo 



( eV 

a 

V^bTk 



■ CT7 



(3) 

In this expression a represents the T = curvature of the 
Kondo resonance with respect to V while 7 describes the 
rate at which the resonance peak broadens and decreases 
in amplitude with increasing temperature. 

Figure [2] shows the values of a and 7 extracted for 
each sample fitted over the range T < Tk/5 and eV < 
fcB'7k/2. We calculate average values of a = 0.051 ± 0.01 
and 7 = 0.107 ± 0.027. Devices with complex 1 had a 
roughly 10% larger average a and 18% larger average 7 
compared to Ceo-based devices. This is due to a few 
samples containing complex 1 which possess exception- 
ally large values of a and 7. These samples may originate 
from conduction occurring closer to the mixed valence 
regime, leading to non-Kondo processes preventing an 
accurate description of the conductance by a universal 
scaling function. Beyond this minor difference, the ex- 
tracted coefficients are essentially molecule- independent. 



Transport measurements (G(T, V) — Gb vs. V) in de- 
vices A, C, and D from Fig. [1] are shown in Figs. [3^ 
[Sb, andlSh, respectively. In Fig. [3ji we plot scaled con- 
ductance, (l-G(T,l/)/G(T,0))/ay, versus (eV/kBTK)'^ 
for the data in these figures, along with the equivalent 
set of traces from device B, for (ey/fcBTk)^ < 2 and 
all available temperatures in the range T/Tk ^ 1.0, us- 
ing the average extracted values of a and 7. We define 
&v = CTa/{l + CT - 1) (T/Tk)2). This plot provides 
a useful qualitative tool indicating the degree to which 
transport in the Kondo regime conforms to the predicted 
scaling function. Data from both devices A and B ad- 
here to the scaling form for (ey/ZcBTk)^ S 0.6 and for 
T/Tk ^ 0.75, above which deviations exceed the con- 
fidence levels. For devices C/D the data conform to 
the scaling function up to (ey/fceTk)^ ^ 0.5/0.55 and 
T/Tk < 0.55/0.45. When data for each device is fit to 
find that device's particular optimal values of a and 7, 
the goodness of fit is essentially the same for every de- 
vice. If we assume the premise of universal scaling such 
that a single (averaged) a and 7 should describe all the 
devices, data from every device with extracted values of 
a and 7 that lie within the average error bars collapse 
onto the scaling curve (black line in Fig. [3ji) up to at 
least (eV/kBTK)^ < 0.5 and for T/Tk < 0.4. 

While the conductance in these devices follows the 
expccted[ll, [H] scaling form, we consider the purported 
universality of the coefficients a and 7 by examining the 
level of accord with other experimental data and theo- 
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FIG. 3: (color online). Conductance as a function of V for 
(a),(b) T/Tk < 1.0 and (c) T/Tk < .65, corresponding to 
devices A, C, and D from Figs, [ij; and[l^, respectively. Gb 
has been subtracted off. (d) Scaled conductance, AG/ay, 
versus (eV/fcaTK)^ for devices A, B, C, and D, where AG = 
(1 - G{T,V)/G{T,0)). The data from (a), (b), and (c) are 
plotted here, along with the analogous data set from device B, 
using the average extracted values of the scaling coefficients. 
The solid black line represents the associated universal curve. 
The colorbar in (b) pertains to all 4 figures. 
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retical expectations. The values extracted from selected 
experimental measurements of nonequilibrium transport 
through a GaAs dot (Tk « 0.3 K) with reasonably 
symmetric coupling (2:1) indicate 0.1 ^ a < 0.15 and 
7 w 0.5. [31, which deviate from our extracted values 
beyond the statistical uncertainties in the parameters. 
Theoretical treatments based on the Anderson (ll. 12. 17| 



and Kondo[13|, [18|, |28j models focus on defining a, 



which is predicted to take on universal values in the lim- 
its of symmetric or strongly asymmetric device coupling. 
Predictions using various approaches are in the range 
0-1 ^ < 0.3, depending on the perturbation method 
used. Our smaller extracted values of a imply that the 
Kondo resonances that we observe evolve more slowly 
than expected with bias voltage given the Tk values in- 
ferred from their equilibrium temperature dependencies. 

There are a number of possible explanations for this 
systematic difi^erence in a. The relative asymmetry of 
the source and drain coupling may be a relevant issue. 
The total level broadening F = Fs -I- Fd , where Fs and 
Fd are contributions of the respective tunneling barri- 
ers as established by the overlap between the local mo- 
ment and the conduction electron states of the source 
and drain electrodes, respectively. In a SMT these val- 
ues are exponentially sensitive to the precise molecule- 
electrode configuration and will generally be different for 
every device. We can infer the asymmetry of the coupling 
from the magnitude of the Kondo resonance as T ^ 0, 
predicted to be (2eV/i)(4FsFD)/(Fs + Fd)^[2^ For the 
SMTs studied here we calculate FgiFn ratios between 
5.2:1 and 182:1 with a median ratio of approximately 
17:1. No systematic difference in coefficient values was 
found for SMTs with increasing coupling asymmetry ra- 
tios. There were also no clear correlations found between 
scaling coefficient values (and thus deviations from the 
scaling curve) and values of Tk, Gq, or Gb. 

It is also possible that the molecular character of the 
devices is relevant to the systematic difference in a. 
The essentially identical observations in Ceo and com- 
plex 1 devices suggest that there is likely a single ex- 
planation that does not depend in detail on molecular 
structure. This argues against higher-spin Kondo states 
or orbital degeneracies, since such physics would differ 
significantly between molecule types. Most theoretical 
treatments assume extremely large (or infinite) on-site 
repulsions relative to the single-particle level spacing, a 
relationship that is not necessarily valid in SMTs. Fur- 
thermore, molecular devices have vibrational de gree s of 
freedom, [HOl These modes have been suggested 31 , S^l 
as relevant to the Kondo regime in SMTs, particularly 
with respect to the unusually weak gate dependence of 
the Kondo resonance, In addition to molecule-specific 
intramolecular vibrational modes, SMTs possess rela- 
tively generic center-of-mass vibrational modes (s^l at low 
energies (~ 5 — 10 mcV) that may be relevant across many 
different molecule types. If the temperature dependence 



of the Kondo resonance in SMTs is modified by vibra- 
tional effects, our extraction of Tk values via fitting the 
temperature dependence may lead to a different effective 
a than in the bare Kondo system. 

We have measured conductance through SMTs in the 
Kondo regime and found that nonequilibrium transport 
at low energies is well described by a scaling function of 
the same functional form as that used for a GaAs quan- 
tum dot in the Kondo regime. A quadratic power-law 
in temperature and bias characterized by two scaling co- 
efficients accurately reflects the evolution of the Kondo 
resonance for SMTs made with either Cgo or complex 1. 
The extracted values of the scaling coefficients deviate 
systematically from those in the GaAs case and from the- 
oretical predictions based on idealized models. These de- 
viations highlight the need for further experimental and 
theoretical investigations of molecular devices and for an 
improved understanding of the limitations of universality 
of Kondo physics in realistic nanostructures. 
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